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Neutrino collective excitations are studied in the Standard Model at high temperatures below 
the symmetry breaking scale in the regime T ^ Mw,z{T) ^ gT. Two parameters determine the 
properties of the collective excitations: a mass scale = gT /A which determines the chirally sym- 
metric gaps in the spectrum and A = M^{T)/2m^T. The spectrum consists of left handed negative 
helicity quasiparticles, left handed positive helicity quasiholes and their respective antiparticles. For 
A < Ac = 1.275 • • ■ there are two gapped quasiparticle branches and one gapless and two gapped 
quasihole branches, all but the higher gapped quasiparticle branches terminate at end points. For 
Ac < A < 7r/2 the quasiparticle spectrum features a pitchfork bifurcation and for A > 7r/2 the 
. . . collective modes are gapless quasiparticles with dispersion relation below the light cone for k <^ 

\ approaching the free field limit for k ^ rUi, with a rapid crossover between the soft non-perturbative 

. to the hard perturbative regimes for k ~ m^.The decay of the vector bosons leads to a width of the 

' collective excitations of order g^ which is explicitly obtained in the limits fe = and k ^ m^A. At 

high temperature this damping rate is shown to be competitive with or larger than the coUisional 
damping rate of order Cf- for a wide range of neutrino energy. 
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> ' I. INTRODUCTION 

00 

Neutrinos are emerging as the bridge between particle physics, astrophysics, cosmology and nuclear physics. Phys- 
ical processes involving neutrinos in dense and hot media are important from cosmology0[ S Hi S S la to the 



, astrophysics of compact stars |8|,|9|,ll0|,lll|. More recently, it has been proposed that Icptogenesis in the early Universe 
lO ' could be the main mechanism that explains the origin of the baryon asymmetry 12, 13, 14]. In this scenario the spec- 
\ trum of neutrinos is an important ingredient that determines the size of the corrections in the transport equations 
(~| • for leptogenesis0|. More refined studies of thermal leptogenesis include finite temperature corrections in the fermion 
O i' propagators [Tsj. which were previously found in vector-like theories (QCD and QED) 16, 17] . 

JL : Neutrino propagation in hot and or dense matter has become very important in astrophysics and cosmology after 
Wolfenstein(18>] studied the effective potential and refractive index of neutrinos in matter and Mikheyev and Smirnov 
(~| , recognized the possibility that matter effects result in resonant flavor oscillations^^. The MSW effect is now rec- 
!• ognized to be the solution of the solar neutrino problem. The importance of medium effects in the propagation of 
neutrinos warrants an active program to study the possible cosmological and astrophysical consequences of novel 
dispersion relations and non-equilibrium aspects of neutrinos in hot and dense media. 
^ \ Early studies of neutrino propagation in hot and or dense media focused on the possible modifications of the neutrino 
dispersion relations in the regime of temperatures relevant for stellar evolution or for big bang nucleosynthesis, namely 
up to ^ 1 — 10 Mev. Neutrino dispersion relations and damping rates in this regime of temperature has originally been 
obtained in ref. 20] up to lowest order in the momentum dependence of the weak interactions in the standard model 
at temperatures well below the electroweak scale. Since then these results have been extended by several authors 
to include the corrections from light scalars|22, to study matter effects in the oscillations of neutrinos in the early 
universe 0, |22, 113, 0] and to a medium that includes nucleonSjkptons and neutrinos [l^ as well as the relaxation 
(damping) rate of neutrinos in a hot and dense medium|^, |2^ |2^ |2^. Damping rates and mean free paths of 
neutrinos in models beyond the standard model have been studied in ref. |28| . 

While all of these thorough studies of the aspects of neutrino propagation focus on temperature (and chemical 
potential) scales suitable for primordial nucleosynthesis or stellar evolution, a systematic study of neutrino dispersion 
relations at high temperatures, near the electroweak scale has not (to the best of our knowledge) yet emerged. 

Motivated to a large extent to explore the possibility of thermal leptogenesis in the standard model, as well as to 
study in general possible non-equilibrium aspects of neutrinos that could be relevant to lepto or baryogenesis|2^ the 
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purpose of this article is precisely to bridge this gap and study the propagation of neutrinos at high temperatures, 
near (but below) the electroweak scale in the standard model. 

Since in the standard model baryon number and CP are not conserved it was conjectured that the cosmological 
baryon asymmetry could have been generated at the electroweak phase transition, if it was of first order. However 
a substantial body of work has revealed that for a Higgs mass larger than about 72 Gev there is a smooth crossover 
instead of a sharp transitionjsoj. The current LEP bound for the standard model Higgs mass mn }^ 115 Gev, all 
but rules out the possibility of a strong first order phase transition and suggests a smooth crossover from the broken 
symmetry into the symmetric phase in the standard model. Supersymmetric extensions of the standard model may 
accommodate a strong first order phase transition, but we will focus our study on the standard model in this article. 

The goal of this article: We study the properties of neutrino collective excitations at high temperatures but 
below the symmetry breaking scale in the standard model. The focus is on obtaining the dispersion relations and 
damping rates of collective excitations of neutrinos in the medium by including the full momentum dependence in the 
contribution from the vector bosons. 

If, as strongly suggested by the numerical analysis 30j , the standard model has a smooth crossover between the 
broken and the unbroken symmetry phases, or even a second order phase transition, the expectation value of the 
neutral component of the scalar doublet diminishes monotonically as the symmetry breaking temperature scale is 
approached from below. In this case the mass of the vector bosons along with the (chiral symmetry breaking) masses 
for the leptons and quarks vanish continuously as the symmetry restoration temperature is approached from below. 

Even when at finite temperature the non-abelian vector bosons acquire a Debye (electric) mass mu ~ gT and a 
magnetic mass m™ ^ g^T 31, 33. IS.^] . an important consequence is that there is a large abundance of vector bosons 
in the medium which enhances processes that are neglected at lower temperatures. 

As it will be discussed in detail below some of these processes lead to a damping rate of collective excitations at 
lowest order in the weak coupling. 

The propagation properties of neutrinos in a medium are obtained from an evaluation of the self-energy which takes 
into account the full propagator of the vector bosons up to one loop order in the high temperature limit. The main 
assumptions that are used in our study are the following: 

• We consider standard model neutrinos, namely left handed massless neutrinos with flavor number conservation. 
The study of neutrino mixing and oscillations will be presented elsewhere |34j. 

• We study the high temperature regime below the symmetry breaking electroweak scale. Tew in which the 
masses of the vector bosons Mw,HiT) « T < Tew- The temperature dependence of the vector boson mass 
arises from the temperature dependence of the expectation value of the neutral scalar in the standard model 
scalar doublet. Under the assumption of a smooth crossover, as seems to be indicated by the lattice data, 
or of a second order transition, this expectation value monotonically diminishes as T ^ Tew from below. 
Since leptons and quarks obtain masses mf(T) via Yukawa couplings to the same scalar that gives the mass 
to the vector bosons the ratios remain almost constant, namely m f{T) / Mw,h{T) ^ m/(0)/MvK,_ff(0) ^ 1. We 
neglect possible logarithmic variations of the weak and Yukawa couplings with temperature via the standard 
renormalization group running. Therefore at high temperature we will neglect the fermion masses in our analysis 
and will consider neutrinos and leptons to be massless. 

• We also consider a CP symmetric thermal bath, thereby neglecting contributions which are proportional to the 
baryon asymmetry of the Universe. 

Brief summary of main results: The main results obtained here are the following: 

• We obtain the effective Dirac equations of neutrinos in the medium from a linear response approach. Under 
the assumption that in the standard model there is either a smooth crossover or a second order transition, 
we consider the regime of temperatures near but below the electroweak scale within which the mass of the 
vector bosons is much smaller than the temperature. In this regime we implement the hard thermal loop (HTL) 
approximationp^ ISlL l32j | and obtain the dispersion relations of neutrino collective excitations in the medium. 
It is argued that perturbation theory and the (HTL) approximation are reliable for T ^ Mw,z{T) ^ gT. 

• In this regime, the properties of the neutrino collective modes are qualitative and quantitatively different from 
those of the fermionic excitations in QCD and QED[i1 [H HI Is^ and depend on two parameters: a chirally 
symmetric mass scale = gT/A and the dimensionless combination A = M^(r)/2 m^T. Collective excitations 
are left handed negative helicity quasiparticles and left handed positive helicity quasiholes and their respective 
antiparticles. For A < Ac = 1.275- •• there are two gapped branches of quasiparticle excitations and one 
gapless and two gapped branches of quasihole excitations. The lower quasiparticle branch and all the quasihole 



3 



branches terminate at end point values of the momentum k determined by A while the dispersion relation of the 
upper gapped quasiparticle branch merges with the free field dispersion relation for k ^ TOj^. For Af, < A < 7r/2 
the quasiparticle spectrum features a pitchfork bifurcation and for A > 7r/2 there are no quasihole branches 
available and only a gapless quasiparticle branch remains. Its dispersion relation is below the light cone for 
soft momenta k <^ and approaches the free field limit for k S> rrii, A with a sharp crossover between the 
non-perturbative and perturbative regimes at k ^ m^. 

• Quasiparticles and quasiholes acquire a width as a result of the decay of vector bosons into neutrinos and leptons. 
This is a novel feature which distinguishes these excitations from those of high temperature QCD and QED and 
has a simple kinetic interpretation. We obtain the damping rate for quasiparticle and quasiholes both at rest 
as well as for fast moving excitations with k ^ my A. These widths are of 0{g^) and it is argued that they can 
be larger than the coUisional width of Cf(7^^|2(l.l2M.l27l| at high temperature for momenta k < MwiT). 

• The domain of validity of the one-loop (HTL) approximation as well as the gauge invariance of the results are 
discussed in detail. 

The article is organized as follows. In section ijl} we obtain the effective Dirac equation for neutrinos in the medium 
implementing a real-time approach in linear response. In section (jIII|l we study the one-loop self energy in the (HTL) 
approximation. In section IjlVj) the spectrum of collective excitations and their damping rates are obtained. In this 
section we discuss the regime of validity of the (HTL) approximation and give a simple kinetic argument that shows 
that the decay of vector bosons implies a width for the collective excitations. Section HVI|) is devoted to a discussion 
of the results and their gauge invariance as well as a comparison between the collisional damping rate of 0{Gj^) and 
the one-loop damping rate. Our conclusions and further questions are summarized in section ()VII|) . An appendix 
provides the technical details for the self-energy. 

II. THE EFFECTIVE DIRAC EQUATION IN THE MEDIUM 

Anticipating a study of the real time dynamics of neutrino propagation in a medium, including eventually mixing 
and oscillations |34|. we obtain the real-time effective Dirac equations for neutrinos propagating in a thermal medium. 
While in this article we focus on standard model neutrinos, the generalization to include mixing, sterile, right handed 
or Majorana neutrinos is straightforward. 

The computation of the self-energy is carried out in the unitary gauge. This gauge reveals the correct physical 
degrees of freedomj35| and previous calculations of the neutrino self-energy in covariant gauges (one of which is the 
unitary gauge) have proven the gauge independence of the dispersion relations [25j. 

The part of the standard model Lagrangian density which is relevant to our discussion is the following 

CsM = (t^) Lva + Clr+d'z+ Ccc + Cnc (II. 1) 

where ^ are the free field Lagrangians for the vector bosons in unitary gauge, namely 

C-w = ~\ {^^,W+ - d,W+) (a^M^-" - d'^W-^') + M^W+W-'' (H.2) 

c°z = -i(a^z,-a,z^)(9'^z"^-a''z^) + iAf|z^z^ (ilb) 

and the charged and neutral current interaction Lagrangian densities are given by 

Ccc = -j= [T^a^L laW+ + la^LiyaW^] (H.4) 



In the above expressions L = (I — 7'^)/2 is the projection operator on left handed states, a is a flavor index and 
g^'^ are the vector and axial vector couplings of quarks and leptons, / stand for leptons and / gcnerically for the 
fermion species with neutral current interactions. 
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The effective Dirac equation in the medium in real time is derived in Hnear response by implementing the methods of 
non-equilibrium quantum field theory described in ref . (361] . In this approach an external source induces an expectation 
value for the neutrino field and this expectation value obeys the effective equation of motion in the medium, including 
self-energy contributions^. This approach allows the study of the real-time dynamics as an mtt^aI value vrohlem 
which will be a distinct advantage to study of neutrino oscillations in the medium. 

To implement this approach we introduce an external Grassmann- valued source that couples linearly to the neutrino 
field via the Lagrangian density 

>Cs = Var\a + ria^a (11.6) 

In presence of this source term, the total Lagrangian density is given by £54/ + -Cs- 

The calculation of expectation values (rather than in-out S-matrix elements) requires the generating function of real 
time correlation functions, which involves evolution along forward and backward time branches|37j. A path integral 
representation of this generating functional is given by 

Z[j+,j-]^ J I?$+I?$-e^/(^[*^J+l-^[*"'J^l) (II.7) 

where we have generically denoted by <I>^ all the fields defined along the forward (-I-) and backward (— ) time branches 
to facilitate the notation. The sources are linearly coupled to these fields in order to obtain the real time correlation 
functions as functional derivatives of the generating functional with respect to these sources. While the sources 
are introduced to obtain the real time correlation functions and are set to zero after the calculations, the external 
(Grassman) source 77 is the same along both branches ± and induces an expectation value of the neutrino field which 
is the same along both branches. See references 36, 37, 38, 39] for further discussion. The equation of motion for the 
expectation value of the neutrino field induced by the external source is obtained by shifting the field 

ryt = xa + ^t ; xa^ii^t) ; (*^>=o (11.8) 

The equation of motion for the expectation value is obtained by requesting that the (^^ ) = be fulfilled order by 
order in perturbation theorv|35 139|]. Implementing this method up to one-loop order we find the following equation 
of motion;36j 




(11.9) 



FIG. 1: One loop diagrams contributing to the neutrino self-energy. The tadpole diagrams (neutral currents) with neutrino 
and lepton loops are not considered in the case of a CP symmetric medium. 

Where SJJ^* {x — x' ,t — t') is the retarded self-energy in real time whose one loop contributions are displayed in fig. 
(Q. In this study we will consider only a CP symmetric medium, therefore we neglect the neutral current one- loop 
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tadpole diagrams with a closed neutrino or lepton loop and only the one loop diagrams with the exchange of charged 
or neutral vector bosons is considered, namely the two top diagrams in figure (|T|). 

The one loop contributions from vector boson exchange to the retarded self-energy are computed in the appendix. 
It is convenient to introduce the space-time Fourier transform of the expectation value and the source 



Xaix^t) ^^T.J dLoXaik,co)e''-'e-^'^' ; rjaix,t) = J ^/a(fc, c.)e*-~ e"^" * (11.10) 

k k 

where we have kept the same name for the variables to avoid introducing further notation. 

Since we are considering the high temperature limit with T much larger than the lepton masses (only the lightest 
leptons are considered) we neglect the fermion masses in the computation of the self energy. As it will become clear 
from the discussion below, this is warranted because, in the high temperature limit the loop momentum is of order T. 
In this case the self-energies are proportional to the identity matrix in flavor space. Using the results of the appendix, 
the equation of motion in frequency and momentum for each individual flavor becomes 

7°cj - 7 • fc + SvF(fc,t^) + Sz(fc,t^) Lxa{k,uj) = -ria(k,uj) . (II. 11) 

The one loop contribution to the self-energies T,w,z{k, oj) are given by the first two diagrams in fig.lpj. The details 
of their calculation is presented in the appendix and their final expressions are given by equations (|A.29|I in terms of 
the dispersive forms (|A.26IA.27|) . 

III. THE SELF ENERGY: CHARGED AND NEUTRAL CURRENTS 

From the results of the appendix, in particular equations ljA.29(l and ljA.30|l we obtain the following explicit expres- 
sion for the self-energies: 

V ft \ f dkn ImT,w,z{k,ko) 

i^W,z{k,Uj) = / ; '■ : (III.l) 



Im^wiKko) = ^/ [l-Nf{p) + m{W,)]mP,q)S{ko-p-Wg)+ ^{-p,-q)S{ko+p + W,)] 



[Nf{p) + m{Wg)]mP,-q)Siko-p + Wg)+ Q{-p,q)S{ko+p~W,)] } (III.2) 



where 



^'^(p, q) = p^ + 2q^ { ^M^^^'^^ 1 (ni.3) 



- (|fc-gl,fc-g) ; - ( 9 ) (III.4) 
= ^q^ + M^iT). (111.5) 

(III.6) 

The contribution from neutral currents is obtained from the above expression by the replacement 

' • ' ; M^[T)^M,{T)^^^ (111.7) 



2COS0„ ' COsOn 

in what follows we will use the standard model value sin^ 9^ = 0.23 for numerical analysis, neglecting possible finite 
temperature corrections. 
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The delta functions in the expression for HIII.2P (and similarly for the neutral current contribution) have a clear 
interpretation: those of the form S{p-\- Wq =F k^) correspond to the process of decay of a neutrino (positive or negative 
energy) into an (anti) lepton and a vector boson and the inverse process of recombination. These processes remain at 
zero temperature but they imply a production threshold that can only be satisfied if the neutrino features a gap in 
its spectrum with a value larger than the mass of the vector bosons. The delta functions of the form 5{Wq — p T ^o) 
correspond to the processes of vector boson decay into a neutrino- (anti) lepton pair and its inverse (recombination) 
process. These processes are only available in the medium and their contribution vanishes in the zero temperature 
limit. 

The form of the imaginary part of the self energies suggests that the full self-energies can be written as follows 



{nil 



The corresponding scalar functions (jy^ ^{k,uj)T<7\^ ^{k,uj) are obtained by projection. Their explicit expression is 
given by 



g'^-K f cPq 1 



[1 - Nf{p) + Nb{Wq)] Qo{p, q) [<5(fco -P-Wq)+ Siko +p+ Wq)] 



2 J {2tt)^ AWqP 
[Nfip) + NtiWq)] Qoip, -q) [5{k^ ~P + Wq) + Siko +P~ Wq)] 



(III.9) 



Imal^iko, k) 



g'^n f dPq 1 



[l-N}{p)+Nk{Wq)]'k-Q{p,q) [5{ko'P^Wq)-5{ko+p+Wq)] 

(III. 10) 



2 J (27r)3 4H/,p 
+ [Nf{p) + Nb{Wq)] k • Q{p, -q) [5{ko -p+Wq)^ 5{ko+p- Wq)] 



Using the dispersive representation we find the following expressions 



d^q 1 



i3 AWqP 



[l-Nf{p)+NbiWq)] Qoip,q) 



1 



1 



[Nf{p)+Nb{Wq)]Qo{p,-q) 



p + Wq ~ uj — ie p + Wq + u + ie 
1 



p — Wq — u! — ie p — Wq + uj + ie 



(III.ll) 



(rw{i^,k) 



r 

2 



d^q 1 



(27r)3 AWqP 
+ [Nfip)+m{Wq)]k-Qip,-q) 



[l-Nf{p)+Nk{Wq)]k-Q{p,q) 
1 



1 



1 



p + Wq — UJ — ie p + Wq + uj + ie 
1 



P — Wq — UJ — ie p — Wq + UJ + ie 



(III.12) 



The expressions for the neutral current contributions can be obtained by the simple replacement l|III.7|) and adding 
both contributions defines the following scalar functions 



a"{u), k) = ayy{u), k) + ct|(w, k) ; (t^(w, k) = <tI^,{uj, k) + (j\{uj, k) 



(III.13) 



An important property of the imaginary parts that will be useful in the analysis of the width of the collective modes 
is the following 



Imcr''(-w,fc) = lxna°{uj,k) 
Im(T"'^(— w, fc) = —lm.a'^{uj,k). 



(III.14) 
(III.15) 



7 



A. Hard thermal loops 

In the high temperature hmit T ^ Mw,z{T) the integrand in the expressions above are dominated by loop 
momentum q ^ T . This is simply gleaned from the high powers of momentum that multiply the Fermi-Dirac 
and Bose-Einstein distribution functions. In this limit, the hard thermal loop (HTL) approximation[l7l Isil l3^ is 
warranted. For frequency and momentum of the neutrino excitations fc, u; <C T the (HTL) "counting" q ~ T ^ fc, a; 
leads to the following approximations 



p^q-k-q ; Wq ^ q + — 

2q 



(III.16) 



where M stands generically for Mw.ziT). 

Ignoring the vacuum contribution, a lengthy but straightforward calculation using these approximations leads to 
the following results for the real parts 



RecrJ^/(fc, w) 



where 



12 M2,(T) 2k 



m,, k 



Tfl 



(III.17) 
(III.18) 



(III.19) 



The dimensionless functions /(w, k) and J(w, k) are more compactly expressed by introducing the dimensionless 
variables 



and the functions 



In terms of which 



my 



k = 



A 



M^{T) ^ g fMw{T) 
2 TO,, T i 



LP(t:J, fc. A; z) 
LM (uJ, fc, A; z) 



2 





k + 


A 

z 


U — 


k + 


A 

z 


w + 


k- 


A 

z 


ID — 


k- 


A 

z 



(III.20) 
(III.21) 



(III.22) 
(III.23) 



/(c.,fc,A) = 



2ze~^ — — 

dz — [iPfcJ, fc. A; z) + LM (uJ, k, A; z)] 

1 — e 



(III.24) 



J(u;,fc,A) 



dz 



2ze" 



1-e- 



■LP{uJ,k,A;z) 



■ LM{uJ, k, A; z) 



(III.25) 



Similar expressions are obtained for the neutral current contributions by the replacement ljIII.7|) . Adding together 
the charged and neutral current contributions leads to 
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lTiLO,k,A) 



JT(w,fc, A) 



J(w,fc, A) 



1 



2 cos^ 
1 

2 cos^ 6„ 



I [uj,k, 



J { uj,k, 



A 



(III.26) 
(III.27) 



We note that for Mw — 



/(w,/c,0) = In 



(III.28) 



J(w,fc,0) = 2 - - In 



(III.29) 



which are the standard results of the (HTL) approximation for vector-hke theories without spontaneous symmetry 
breaking (QCD-QED)^!!!!!!!. The following properties of these functions will be important in the discussion of 
the collective modes below, 



/(-w,fc, A) 
J(-w,fc, A) 
lieu, -k, A) 

J{LJ,-k,A) 

^(c^,0,A) 



-/(w,fc,A) 

J{uj,k,A) 

-/(w,fc,A) 

J{uj,k,A) 





(III.30) 
(III.31) 
(III.32) 
(III.33) 
(III.34) 



IV. DISPERSION RELATIONS AND WIDTHS OF COLLECTIVE EXCITATIONS 



The form of the self-energies (jIII.8|l allow to write the homogeneous Dirac equation (setting to zero the external 
Grassmann sources i] ) in the form 



A+ik)D^{u;, k) + A_(k)i?+(c^, k) Lxa{k,Lu) = 



where 



A±(k) = i(7°T7-k)-|^(lTMk)7') 



here /i(k) is the helicity operator, and 



D+{LU,k) = LU - k + (7"{k,ui) - a^{k,uj) 
D_{uj,k) = uj + k + a"{k,uj) + (T'^{k,uj) 



Therefore, the propagator for left handed fields is given by 



(IV.l) 



(IV.2) 



(IV.3) 
(IV.4) 



iS{uj, k) = 



A+(k) , A_(k) 



D+{Lu,k) D^{uj,k) 

The neutrino spectral function is determined from the propagator above, and is given by 



(IV.5) 



p(a;,fc) = --ln\[iS{Lo + ie,k)] = A+(k)p+(w, fc) + A_(k)p_(w, fc) 

TT 



(IV.6) 
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where 



Im D±{lu + ie, k) 



'R.eD±{uj + ie,k) + lYaD±{oj + ie,k) 



(IV.7) 



The poles of the propagator along the real axis in the complex lj plane correspond to physical stable excitations, 
while complex poles very near the real axis in a second (or higher) Riemann sheet in this plane, describes quasiparticles 
or resonances. Propagating quasiparticles are characterized by narrow width resonances, namely their widths must 
be much smaller than the real part of the (complex) pole, so the decay rate of the quasiparticle is much smaller than 
its oscillation frequency. 

The position of the (quasi) particle poles are Lo±{k) which are determined by the conditions 



ReD+{uj+{k),k) = 
ReD.{Lj.{k),k) = 



(IV.8) 
(IV.9) 



// these functions vanish linearly near their zeroes then near the resonances the spectral density can be approximated 
by a Breit-Wigner form, namely 



p±{ui,k) 



Z±{k) 



r±(fc) 



(c.-c.±(fc)) +rl{k) 

where the residues Z±{k) and widths (damping rates) T±{k) are given by 



(IV. 10) 



z±\k) 



dD±{uj,k) 



duj 



(IV. 11) 



^±(fc) 



r± (k) = Z±{k) Im D± {uj± (k) , k) 



(IV.12) 



The spinor wave functions of the resonances (quasiparticles) are obtained from the homogeneous effective Dirac 
equation in the medium (jlV.ip by considering only the real part of D±{uj, k), since the vanishing of the real parts 
defines the quasiparticles. Because neutrinos are left handed fields the eigenspinors corresponding to the solutions of 

D-i-{uj+{k), k) = must obey 



il + hik))Lx{Lo+ik),k)^0 



(IV. 13) 



Namely the spinor eigenstate is a left handed, negative helicity state, just as the usual neutrino field. We will refer 
to the positive energy solutions of Dj^{LU^{k), k) = as quasiparticles. Similarly, the spinors solutions corresponding 
to the roots of D^{uj^{k), k) = must obey 



il^h{k))Lx{iU-{k),k) = 



(IV. 14) 



Corresponding to a left handed neutrino state with positive helicity. Following ref.|l7|. we will refer to the positive 
energy solutions of £'_(w_(fc), /c) = as quasiholes. These states feature the opposite chirality-helicity assignment of 
the vacuum neutrino states. 

As it will be discussed in detail in section llV)l . the regime of validity of the perturbative expansion is restricted 
to Mw,z{T) ^ gT in which case m^/Mw,z{T) ^ 1. Therefore the first terms in equations ljIII.17IIII.18|l provide a 
perturbative correction to the coefficients of uj and k in D±(uj, k). Since the focus is to explore the non-perturbative 
aspects of "soft" collective excitations we will ignore these terms in what follows. The final expressions for D±{uj,k) 
are given by 



D^{u!, k) — to ■ 



2k 



f- [/t(w, k) + Jt(w, k)] + i [/mcr"(w, k) - Ima\uj, k)] 



(IV. 15) 
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(IV. 16) 



The real and imaginary parts of D±{lo, k) feature the following symmetries 

ReD+{-uj,k) ^ -ReD^{uj,k) 
ImD-^-{—uj, k) — ImD^ (uj, k) 



(IV. 17) 
(IV.18) 



The dispersion relations of quasiparticles and quasiholes are obtained from roots of the real part of the self-energy, 
namely 



ReD±{uj±ik),k) = 



(IV.19) 



The symmetry relation (|IV.17p indicates that we only need to find the positive roots in eqn. (|IV.19|) . therefore the 
roots of ReD±{uj, k) — will be the pairs {uj+{k), —uj-{k)); (a;_(fc); —uj+{k)) respectively, with uj±{k) > 0. 

Antiquasiparticles and antiquasiholes have negative energy and quantum numbers opposite to those of quasiparticles 
and quasiholes. The interpretation of the excitation spectrum is as follows [l7j: the positive energy roots of -D+(ti^, k) 
correspond to negative helicity left handed quasiparticle states of energy w+(fc) and spectral weight Z+{k) while the 
negative energy roots are left handed anti-quasihole states of negative helicity of energy —u}^{k) and residue Z^{k). 
The positive energy roots of D_(uj,k) correspond to left handed positive helicity quasiholes of energy Lu_{k) and 
spectral weigth Z_(k) and the negative energy roots correspond to left handed positive helicity anti- quasiparticles of 
energy — w+(fc) and residue Z+{k). The relation between the imaginary parts, eqn. I|IV.18|) states that quasiparticles 
and antiquasiparticles have the same width and so do quasiholes and antiquasiholes, this is obviously a consequence 
of the underlying CPT symmetry. 

The table IdcIow summarizes the properties of the collective excitations 



Identity 


E 


h 


Z 


r 


ik) 


quasiparticle 


uj+{k) 


-1 


z. 


-ik) 


r_ 


h(fc) 


antiquasiparticle 


-OJ+{k) 


1 


z. 


-ik) 


r_ 


h(fc) 


quasiholc 


uo^{k) 


1 


z_ 


-ik) 




-ik) 


antiquasihole 


-UJ^{k) 


-1 


z 


-ik) 


r_ 


ik) 



A. Gaps in the spectrum of collective modes 

The position of the gaps in the spectrum of collective excitations are determined by the roots of the equation 



ReL»±(w,fc = 0, A) = 



(IV.20) 



where we have used the property J{uJ, k = 0) = (see eqn. (IIIL34|l ') which implies that D^iuo, fc = 0, A) = D^{llj, k 
0, A). General features of the gap equation l|IV.20|l can be understood from the expression 



D+{uj,0,A) = m^ i ^ + 2 



, \ ' cos^ 6 



(IV.21) 



where 



H{u;,A) = 




(IV.22) 



and V stands for the principal part. 
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For A = 0, H{uj,0) = ~-2V{l/uj) which leads to the solution of the gap equation 

.(fc = 0)^±^.(l + ^)^ (IV.23) 

We note that because V{l/uj) is an odd function of uj vanishing at a; = (by definition of the principal part) there 
are three roots of the eqn. _D+(w,0, 0) = 0: uj — and the two roots corresponding to the value of the gap (|1V.23II . 
The gapless excitations associated with the root oj = will be explored in subsection IjIV C|) below for general values 
of A including A = 0. 

For A 7^ the function D+{uj = 0, k — 0, A) — because it is an odd function of uj. Furthermore, for A ^ 1 and 
u; w 0, it follows that 

H{co,A)^^ (IV.24) 

therefore, for A ^ 0, D+{uj, fc = 0, A) vanishes at w = 0, is positive for u; > and rises sharply from the origin with a 
slope oc 1/A^. For small A the function H{uj, A) eventually becomes negative for large enough oj. This is because the 
integrand is sharply peaked at z ~ 2, therefore, for sufficiently small A and sufficiently large lu the region in which 
the denominator is negative dominates. Since for large uj, D^{ijj, A) « there must be yet another sign change and 
the function D+{uj^A) must have at least three zeroes^ namely three roots. One of the roots corresponds to oj — Q 
which determines a gapless branch of collective excitations, however for small A one of the roots must be at ZU ^ A 
while another root should be at (1 + 2 cos^^ e ) ^ • Therefore for A 1 there are three branches, one emerging 
from the origin is gapless and two other branches that feature a gap in their dispersion relations. 

However for larger A the region in which the denominator becomes negative necessarily corresponds to z >> 2 and 
the integrand is strongly suppressed in which case the function H(u, A) will be always positive and no root will be 
available. This behavior is clearly displayed in fig. Q where we used the standard model value sin^ 9^ — 0.23. 



8.0 




-6.0 I — ^ ^ ' ^ ' ^ ' ^ 1 ^ ' 1 

0.00 0.50 1.00 1.50 2.00 2.50 3.00 



FIG. 2: Roots of D{uj, fc = 0, A) for A = 0.2, 1, 1.3, 5 respectively. The intersection of these curves with the horizontal axis 
determines the values of the gap in the spectra of collective modes. 

This figure shows that there is a critical value Ac such that for A < Ac there are two non-vanishing roots which 
coalesce at A = Ac while for A > Ac the only root available is w = 0. We find numerically that the critical value is 
given by 

Ac- 1.275- •• (IV.25) 

In order to resolve whether the gapless branch describes quasiparticles or quasiholes, we must find the k dependence 
of the gapless branch. 
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A straightforward but lengthy calculation shows that to linear order in k 



D±{u:, k, A) = D+(c., 0, A) T fc 1 + S{uJ, A) + S{iJ, 



with D+{uj, 0, A) given by eqn. HIV.2ip and 



S{lJ,A) = -— / dz 



Stt^ Jq 1 — e 



2 cos^ 6w cos^ 



Therefore the dispersion relations for small k near the gaps with value uj±(0) are given by 



(IV.26) 



(IV.27) 



with the group velocities c± given by 



uj±{k) ~ w±(0) + c±k + 0{k^) 



(IV.28) 



± 



i?'(w±(0),0,A) 



l + S{ZJ±{Q),A) + 



2 cos^ 



A 



cos^ 9„ 



(IV.29) 



Therefore for sufficiently small k all branches of collective modes behave linearly with fc. The study of the group 
velocities for each branch necessarily has to be done numerically, however of particular interest is the group velocity for 
long- wavelength excitations of the gapless branch. A straightforward calculation for uj±{0) = leads to the following 
result 



c± =± 



A2+3£i 



(IV.30) 



where ± refer to the roots of Z?±(a;,fc, A) respectively. As discussed above, quasiparticle and quasiholes are positive 
energy roots of D±{lo, fc, A) respectively. Hence, the expression for the group velocity for long wavelength excitations 
on the gapless branch (jIV.30p indicates that for A < 7r/2 positive energy excitations on the gapless branch are 
quasiholes^ while for A > 7r/2 gapless collective modes are quasiparticles. Therefore the gapless branch changes 
identity at A = 7r/2. 

In summary, the spectrum of collective excitations always features a gapless branch for any value of A. For A < 7r/2 
these gapless collective modes are quasiholes, and for A > 7r/2 they are quasiparticles. 

For generic A < Ac ~ 1.275- • • the quasiparticle spectrum will feature at least two gapped branches, while the 
quasihole spectrum will feature two gapped and one gapless branch. For A ^ 1 the lowest gapped branch begins at 
|w(0)| ~ A while the highest one begins at 10^(0)1 ^ 'm„{l + 27ES^^)^ ■ ^ > ""/^ quasiparticle spectrum only 
features a gapless branch and no branches of collective excitations remain for quasiholes. 

The dispersion relations are linear for small momentum fc in each branch. However, the gapless branch corresponds 
to roots of D_{uj, fc. A) for A < 7r/2 while it describes the only available, gapless branch of D^{uj, k, A) for A > 7r/2. 
A detailed numerical analysis of the quasiparticle and quasihole collective modes is presented below. 



B. Quasiparticle spectrum 

The dispersion relation of quasiparticles corresponds to the positive roots of Re D+{uj+{k), fc. A). The study of the 
dispersion relation in the full range of fc and A must necessarily be carried out numerically. However, before doing so, 
we can gain insight into the nature of the dispersion relations by considering the simplified case of A = 0. While this 
case is outside the domain of validity of the perturbative expansion because of the restriction that Mw,ziT) ^ rui, 
(see section (©), the study of this case highlights several relevant aspects of the dispersion relation and also serves 
as comparison to the more familiar results in the literature in the case of unbroken gauge theories in the (HTL) 
approximation ^3 1^ . 
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In the case A = the function D+iuj, fc. 0) is famihar from the (HTL) approximation for fermionic excitations in 
unbroken vector-Hke gauge theories[T3, Is^, |33| 



C 

D+{lu, k, 0) = m^< U — k -\ — = 

2fc 





In 


w4 


k 










+ 2 


Lv k'' 




UJ — 


k 





C = 1 



1 



2 cos^ 0,, 



(IV.31) 



This function is negative for w ~ but becomes positive for uj ^ k and features only one positive root for any value 
of k. The gap is given by 



2 cos^ 6*,, 



(IV.32) 



The group velocity for long-wavelength excitations around this gap is obtained from the general expression H1V.29() 
by setting A = and the value of the gap (|IV.32|I . A straightforward analysis yields c+ — 1/3. 

This situation is very different for A 7^ 0. The discussion on the values of the gaps in the spectrum of collective 
excitations in the previous subsection (|IV A|) as well as the dispersion relation of long wavelength excitations in the 
gapless branch indicates that there are three distinct regions of the parameter A with qualitatively different behavior. 



A < Ac 



1.275- 



In this region the quasiparticle spectrum features two gapped branches and no gapless branch. We denote the 
dispersion relation for the lowest branch by uj^{k) and and that for the higher branch w>(fc) respectively. 

For A ^ 1 the gap for the lowest branch is w<(0) ~ A whereas the gap for the highest branch is w>(0) ~ 

(see fig. For increasing values of k the roots corresponding to the largest gap move to larger 



values of w, while the roots corresponding to the lowest gap move to smaller values of w. The spectrum for the lowest 
branch terminates at an end point k^ which numerically is found to be k^ < A/2. For k > k^ there are no further 
roots corresponding to the lowest branch and only the highest branch is available. For k 3> rrii, the dispersion relation 
along this remaining branch is found to be 



UJ+{k)^k+-^[l 



2 cos^ 



(IV.33) 



The positive energy roots for A ^ 1 are displayed in fig. ||2J). 



o 

II 

< 

3 




FIG. 3: Roots of D+{iu, k,A = 0.2)/m„ for = 0.03, 0.1, 0.3, 0.5, 1 respectively. 
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FIG. 4: Roots of D+{uj,k,A = l)/m^ for = 0.2,0.3,0.4,0.5 respectively. 



For larger values of A < Ac the behavior is qualitatively similar but with some quantitative differences. The 
positive energy roots of D-^-{u!, k,A — 1) for several values of k are displayed in fig. 

This figure clearly reveals that the branch with the largest gap remains for any value of fc, whereas that corresponding 
to the smallest gap terminates at an end point value of momentum which numerically is found to be roughly fee < A/2. 
For k > ke there are no real roots for the lowest branch. 

The weight of the (weakly damped) collective modes to the spectral function is determined by the residue Z^{k, A) 
given by eqn. H1V.11|I evaluated at the position of the root, ijj_|_(fc,A), which determines the dispersion relation of 
collective excitations on a particular branch. 

The analysis in the previous subsection shows that the function D_|_(a>, fc, A) rises sharply for small w with a slope 
~ 1/A^ therefore for A ^ 1 the weight of these collective modes with dispersion relation corresponding to the root 
A) is 

Z+(fc,A)-A2<l (IV.34) 

Hence for A ^ 1 the branch of collective modes with dispersion relation a;^(fc) has negligible spectral weight. 

We find that for k 3> m^, the asymptotic behavior of the collective modes in the upper branch is similar to the 
(HTL) limit of unbroken vector- like theories |l7l l3ll l33| . For fast moving quasiparticles corresponding to the upper 
branch in the limit k ^ the dispersion relation is given by eqn. ljIV.33p and the spectral weight is found to be 

A) - 1 - §(l + ^) 1" ) + ■ ^ ^ (IV.35) 

The dispersion relations for both branches together and for the lowest branch separately are displayed in fig. (jSJ 
for A = 1. 



2. Ac < A < 7r/2; a pitchfork bifurcation of the spectrum. 

As discussed above, for A > Ac = 1.275 • • • the only solution of the gap equation is w = 0, however for A < 7r/2 
there is no gapless positive energy solution with linear dispersion relation for small k quasiparticles. Therefore the 
following question arises, what is the nature of the quasiparticle spectrum for Ac < A < tt/27. A numerical study 
reveals a remarkable answer: for a given value of A in this region there is critical value of fc, denoted by fcc(A) for 
which two roots of the equation 13+ (w, fc, A) = emerge with a pitchfork bifurcation in the spectrum. This result is 
depicted in fig. © for A = 1.4. 
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FIG. 6: Pitchfork bifurcation of the quasiparticle spectrum for A > Ac. The function D+{u!, k, A = 1.4)/mi, vs. co/m,^ for 
k/m^ = 0,0.1,0.2,0.3,0.5. 



Note in the figure ^ that there are no roots for D+{uj,0, A ~ 1.4), however two roots emerge continuously for 
fc > fcc ~ 0.1 m^. 

For A — Ac the critical value of the momentum vanishes and becomes non-zero continuously for 7r/2 > A > Ac. 
The critical value fee (A) determines the origin of the pitchfork bifurcation. One of the branches of the pitchfork 
bifurcation moves towards larger values of the frequency but the other towards smaller values, eventually terminating 
at an end point which for the case of the figure ^ is kg ^ 0.5 m,^. For large values of the momenta we find that the 
roots determine the same dispersion relation as for the previous cases, which is given by eqn. ljIV.33|) and the spectral 
weight is also given by eqn. (jlV.35|) . 

While this phenomenon of a pitchfork bifurcation in the spectrum is a noteworthy aspect of neutrino collective 
excitations in the medium, the narrow window in the parameter A within which this phenomenon emerges is rather 
restricted and without a particular significance within the standard model. 



3. A>Tv/2. 



For A > tt/2 only the gapless branch of D^{llj, k, A) is available and the bifurcated spectrum ends. This is depicted 
in fig. Q which displays the roots near the origin for small values of k/m^. These roots move continuously towards 
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larger values of the frequency for larger momenta. 
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FIG. 7: Emergence of the gapless branch for D+(Lj,fc, A) for A > 7r/2. 



The dispersion relation for A > Ac is displayed in fig. Q for A = 3. For fc ^ TO,yA the dispersion relation is linear 
but below the light cone and the group velocity agrees with the result (|1V.30|I . whereas for k ^ m^A the dispersion 
relation is found to be 



(IV.36) 



which is the same as for the large momentum limit of the previous cases. Thus the spectrum interpolates continuously 
between the soft non-perturbative region k <C rrii, and the hard perturbative region k 2> with a rapid crossover 
at fc w m,,. 




1 



FIG. 8: Quasiparticle dispersion relation for A = 3. 



The large momentum limit of the spectral weight in this case is also given by eqn. ljIV.35|l . 
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C. Quasi- hole spectrum. 

Just as in the previous section, we begin the analysis by revisiting the case A = 0. In this case 



c 






In 


uj + k 












- 2 


2I 








uj — k 





C = l 



1 



2 cos^ 9„ 



(IV.37) 



For UJ = this function is positive and for tU ^ fc it has the asymptotic behavior D-{LU,k,0) ^ lu. However for 
Lo ^ k the logarithm gives rise to a sharp downward spike and the function becomes negative near the light cone for 
any value of k. Therefore the function features two roots one above and one below the light cone. The dispersion 
relation for small k is obtained from cqn. (|IV.26|I for A = 0. 

As fc ^ the root above the light cone approaches the gap ljIV.32p from below with group velocity c_ — —1/3 
and the one below the light cone approaches the origin from above with group velocity c_ = 1/3. This is precisely 
the gaplcss branch, which is the solution of the gap equation for any A, including A = as mentioned in subsection 



This gapless branch has not been discussed in previous studies of the (HTL) approximation in unbroken gauge 
theories_17, 31, 32] because it is below the light cone and is strongly Landau damped in those cases. We will discuss 
below the width of collective excitations in this branch (see section (|IVD|l ). 

This simple analysis combined with the linear dispersion relation for long wavelength collective modes in the gapless 
branch with positive energy roots of Z)_(a;,fc,A) indicates that the gapless branch is a solution of Z)_(ti;, fc. A) — 
for A < 7r/2. Positive energy solutions along this branch correspond to gapless quasiholes. For A <C 1 the function 
D_(a;, fc. A) has a slope (x 1/A^ for both the gapless branch and the branch with the smallest gap ~ A. Therefore 
the spectral weight for these branches is of order A^ ^ 1 and therefore negligible with respect to that of the branch 
with gap ~ TOi.(l + 2col^eJ ^- 

For A = the origin of the two roots, above and below the light cone, is traced to the logarithmic singularity at 
ijj = k \ii eqn. (|IV.37|) . However inspection of the functions LP and LM given by eqns. (|III.22IIII.23|I reveals that 
for any A ^ this singularity is screened. For small A a similar situation featuring two roots is expected for small k. 
However for larger A and k the screening of the light cone singularity will not allow the function to become negative. 
This feature is clearly displayed in fig. ^ for A = 0.2. 
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FIG. 9; Roots of D-(uj,k,A = 0.2)/m^ for = 0.1,0.2,0.3,0.5, 1.0, 1.5 respectively. 



Inspection of this figure reveals the presence of three branches as anticipated, two gapped branches and one gapless 
branch, which is manifest as the root closest to the origin for small values of k in the figure. The gapless branch 
terminates at an end point value of the momentum ke- Fig. (jSJ shows that for larger values of k there are no available 
roots (see for example the curve for k/rrii, = 1.5 in the figure). Therefore for larger values of the momentum all 
the branches terminate at different end point values that depend on A. It is found numerically that this behavior is 
consistent for all values of A < A^: one gapless and two gapped branches, each branch terminates at a different end 
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point value of the momentum which depends on A. In particular the end point value for the gapless branch diminishes 
continuously as A Tr/2. Figure (|10|1 displays the function D_{uj, fc, A = 1) as a function of oj for several values of 
k. Three branches are clearly displayed in this figure as well as the termination of these branches at particular values 
of the momentum fc, for example only the gapless branch remains for fc > O.Sm^ and finally no roots are available for 
k > OAnii,. The quasihole spectrum features end points for all the branches and there are no propagating quasihole 
states for fc ^ m^. 



2.00 




0.00 0.50 1.00 1.50 2.00 

co/niv 



FIG. 10: Roots of D„(a;,fc,A = l)/m^ for = 0.2,0.3,0.4,0.5 respectively. 



Figure (|11|) display two relevant cases. The left panel shows the function D_((jj,fc,A) for Ac < A < 7r/2. This 
figure shows that the gapped branches have disappeared and there is a rather small window of (small) momentum 
within which there are gapless quasihole states. Finally the figure on the right corresponds to A > 7r/2, which clearly 
reveals that there are no available roots and all of the quasihole branches have disappeared, including the gapless 
branch, consistently with the previous discussion. 




0.00 0.50 1.00 1.50 2.00 I 

0)/mv 



0}/mv 



FIG. 11: Roots oi D-{u},k,A = 1.4) /m^ and D-{uj,k,A = 1.8)/m^ for = 0.1,0.2,0.3,0.5,1,1.3 respectively. The arrows 
indicate increasing values of k/m^ 



In summary, the spectrum of quasihole states features three branches for A < Ac, one gapless and two gapped 
branches, each branch terminates at particular end-point values of the momentum which are functions of A. For 
Ac < A < 7r/2 the gapped branches disappear and only the gapless branch remains but terminating at an end point 
value of the momentum which is rather small. For A > 7r/2 there are no quasihole states. Only for A ^ 1 there 
are quasihole states for fc ^ m^, this is consistent with the existence of quasihole roots for arbitrarily large fc in the 
A = casein m. 
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For generic values of A, quasihole states are non-perturbative, soft collective excitations available only for momenta 
k < m^. 

The spectral weight of the quasihole excitations is related to the slope of the function £)_ {uj, k, A) as a function of ui 
by cqn. lIV.lip . In the general case the spectral weight must be found numerically, however, inspection of the figures 
in this section clearly reveals that the spectral weight associated with the gapless excitations is much smaller than 
that of the collective modes on the gapped branches. For oj <C rrii, the slope of the fmiction D-{uj, k, A) is « 1/A^ 
for A <C 1. Therefore gapless collective modes yield smaller contributions to the spectral density than the gapped 
excitations. 

D. Damping rates of collective excitations: 

The damping rates or widths of the quasiparticle and quasihole excitations are given by eqn. (|IV.12(I where Z±{k) 
are the residues given by eqn. ljIV.ll|) , and the imaginary parts are given by 

lmD±{uj±{k),k) = [lina'^{k,uj)Tlmali,{k,u)] + [lm(T|(fc, tj) T Imcri(^, w)] (IV.38) 

The imaginary parts Ima^ ^{k, uj); Imcrjy ^(fc, oj) are given by the expressions ljIII.9pllII.10|) for the charged current 
interactions and similar expressions obtained via the replacement (IIII.7|) for the neutral current interactions. An 
analysis of the delta functions in these expressions determines the region of support of the imaginary parts and 
thereby establish whether there is a non-trivial width for the collective excitations with dispersion relations uj±{k). 

• Region of support for S{uj — p — Wq): this delta function corresponds to the "decay" process ly W + 1 
for the charged current and v Z + V for the neutral current. Clearly the neutral current contribution 
cannot be fulfilled. For the charged current contribution, the quasiparticle and quasihole could decay provided 
nil, > Mw{T). In principle the kinematics for this process could be satisfied since near the "transition" the 
mass of the vector boson may be small. However this would entail that 

gT > Mw{T) (IV. 39) 



which, however, contradicts the bounds (|V.4|1 . (|V.5p which determine the domain of validity of the perturbative 
expansion, the latter one being the most conservative. Therefore within the domain of reliability of the pertur- 
bative expansion invoked in this study, the kinematics resulting from this delta function cannot be fulfilled on 
the quasiparticle or quasihole mass shell. The delta function 6{u! + Wq + p) has support only for uj < and 
corresponds to the decay of a negative energy neutrino. 

• Region of support for (5((jJ +p~ Wq): this delta function corresponds to the decay process W ^ v + 1 (the 
term with 5{uj — p + Wq) corresponds to the decay into antineutrino-lepton) and it has support on the neutrino 
mass shell for Mw{T) > rrii, for a massless lepton, or Z v + i/ for Mz{T) > 2m^. Either of the bounds (|V.4p . 
IjV.Sp which determine the domain of validity of the perturbative expansion guarantees that these processes are 
kinematically allowed in the region of validity. Therefore this analysis leads to the conclusion that the decay 
of the vector bosons into neutrino-lepton pairs (charged current) or neutrino pairs (neutral current) lead to a 
width for the neutrino collective excitations. These processes are depicted in fig. H12|l . 

The fact that the decay of vector bosons into neutrinos leads to a width for the neutrino states can be easily 
understood from a simple kinetic argument. Consider the kinetic equation for the time evolution of the neutrino 
distribution function. This equation is obtained by a gain minus loss argument as follows. Consider first the charged 
current interaction: W-decay, W ^ vl, leads to a gain in the neutrino distribution and the recombination v + l ^ W 
leads to a loss. The gain and loss terms are obtained from the transition probability per unit time for each process 
leading to the following kinetic equation 



dn^ik) f d^q \Mw{q)\^' 



dt J (27r)3 AWqp 



(1 - n,{k)){l - 7ii{p))Nw{Wq) ~ n,{k) ni{p) (l + NwiWq)) 



S{Wq-p~uj{k)) (IV.40) 



where lAlwl'?)!^ is the transition probability matrix element, p = k — q and n; , N\y are the distribution functions 
for leptons and for vector bosons respectively. A similar contribution arises for the neutral currents which can be 
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FIG. 12: Decay of vector bosons. W ^ v + I via the charged current or Z ^ vu via the neutral currents. The blob in the 
external neutrino line corresponds to a quasiparticle or quasihole. 



simply obtained from the expression above by the corresponding replacements and ni —^ n^. The hnearized kinetic 
equation around the equilibrium distributions, is obtained by writing n^ik) = Nf{k) + S ni,{k) with Nf(k) the Fermi- 
Dirac distribution function and keeping the lepton and vector boson distribution functions to be those of equilibrium, 
namely ni{p) = Nf{p) ; N\Y{Wq) = Ni,{Wq). Keeping only linear terms in S n^{k) one obtains the kinetic equation 
in the relaxation time approximation 



d 5n^{k) 



-Sn^{k) 



d\ \Mw{q)\' 



dt ^ ' J (27r)3 AWqP 

Thus we see that the decay rate of the distribution function 



Nk{Wq)+Nf{p) 



5{Wq - p - uj{k)) . 



(IV.41) 



d^q \MwiqW 



(27r)3 iWqP 



Ni,{Wg) + Nf{p) 



S{Wq~-p-uoik)) 



(IV.42) 



is identified with last terms in the imaginary parts ljIII.9|)llTl.lO|) . The damping rate for quasiparticles T{k) is simply 
related to the relaxation rate of the distribution function as 7(fc) = 2r(fc) JJ,, J6., .32^ ^] . The fact that the decay 
of a heavier particle implies a width for the light collective modes in a thermal bath has already been recognized in 

refs.mii^. 

A simple analysis of the region of support of d{Wq — p — uj) in ljIII.9|)ITTl.lO|) shows that for oj above the light cone 
(lu > k) the loop momentum integral is restricted to the range < q < where 



M2,(r)-(u; + fc)2 



2{uj + k) 



M^{T)-{u-kf 
2{lo - k) 



(IV.43) 



We find that the charged current contribution to the imaginary parts are given by 



ImD^ {uj+{k),k) 
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k-Q{p,~q)) Nf{p) + miWg) 



(IV.44) 



with 



p = Wq~uj+{k) ; q-k 



and 



q^ — k"^ — p^ 
2qk 



(IV.45) 



2 /.9+(w_(fc)) 1 _ 

lu,D^{Lo.{k),k) = -f- / l^(Qo(p, -9) + k • Q{p, -q)) Nfip) + N,{Wq) (IV.46) 
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with 



P - 



W,~LU-{k) ; q-k: 



2qk 



(IV.47) 



The contribution from neutral current interactions is found from the expressions above by the replacement Ijlll. 
The damping rate can be obtained numerically in the general case, but analytic progress can be made in two limits: 
quasiparticle (quasihole) at rest, namely fc = 0, and for fast moving quasiparticles T 3> fc 3> m^. As discussed in the 
previous section, only for A <g; 1 are there quasihole states for k ^ m^. 

Furthermore within the domain of validity of the perturbative expansion determined the tighter bound ljV.5|) it 
follows that Mw{T)/m,y ^ 1 which results in several simplifications. 

Quasiparticle and quasihole at rest: for fc = the damping rates of quasiparticle and quasiholes are the same. 



It can be obtained directly from the expressions for the imaginary parts ljIII.9(lllTl.lO(l by recognizing that for fc = 
the contribution from (|III.10|I must vanish by rotational invariance or alternatively by computing the integral in the 
expressions above in the limit ^ . The final result for the damping rate is 
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(IV.48) 



with 



F[M] = 



^2 



[Nf{q*)+N,{W{q*))] 



M2 -ul ujI 
^ ; W{q*)^ ^ 



2cj„ 



2uJn 



(IV.49) 
(IV.50) 



where Wg is a non-zero solution of the gap equation (|1V.20() and Z[ujg\ is the residue corresponding to zero momentum 
quasiparticles or quasiholes. 

Since Ug < and in the domain of validity of the results Mn/,z(T) 3> m^, the expression for the function F[M] 
m eqn. (jIV.49|) simphfies to 



F[M] 



M 

ijJn 



sinh 



2u„t\ 



(IV.51) 



Therefore the width for quasiparticles and quasiholes at rest is given by 



2 uriij 



1 



sinh 



2cos6 0^ sinh[^ 



(IV.52) 



For A ^ 1 the gap for the lowest gapped branch is a'< ~ rrii, A and for the upper branch is a;> ~ rrii, while the 
residues are Z[u)g\ ~ A^ for the lower gap and Z[ujg\ ^1/2 for the upper branch. 
Hence it follows that for A ^ 1 



r(w<) 



A ujy 



(IV.53) 
(IV. 54) 



Therefore for A <C 1 the lowest branch is strongly damped while the upper branch is weakly damped. 
For A '--^ 1 ; A < Ac both the lower and upper branch have gaps of order (up to factors of order one, see for 
example fig. and Z[ijjg\ ^ Hence for A < Ac but A ~ 1, 



(IV.55) 
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and collective excitations at rest on either branch are strongly damped. 

Fast moving quasiparticles : The other limit that is relevant for a comparison with the perturbative results and 
tractable analytically is that of A > 1 and fast moving quasiparticles with M{T) ^ k ^ A. In this limit only 
quasiparticles are available since the spectrum of quasiholes terminates at an end point of order fcg ^ m^A for not too 
small A (see figs. H9lll|l '). For k ^ rrii^A the quasiparticle dispersion relation is near the light cone and given by cqn. 
pV.33p and the spectrum of quasiparticles is perturbatively close to the free field spectrum, therefore a comparison 
to the perturbative results is meaningful. 

In this limit q^/T ^ 1 and q- ^ M{T) the upper limit of the integral can be taken to (7+ ^ oo. The remaining 
integrals are now straightforward and after combining the results from charged and neutral currents we find 



4fc 



1 + 
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2 cos4 e„ 
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Thus in the limit k 3> mriuA quasiparticles are weakly damped since 



r+(fc)/fc^A^ln 



< 1 



(IV.57) 



V. DOMAIN OF VALIDITY OF THE HTL APPROXIMATION: POSSIBLE CAVEATS. 

In the (HTL) approximation we have used the bare propagators for the internal fermion and vector boson lines. 
Since the (HTL) approximation relies on loop momentum g ~ T the region of validity of the (HTL) approximation 
is determined by the region in which the self-energy corrections to the internal lines are perturbative with respect 
to the loop momentum scale. The study of the previous section indicates that vector boson exchange will lead to a 
radiative correction to the internal fermionic line yielding a mass scale rrii, ~ gT . Thus the self energy correction to 
the lepton or neutrino internal lines is indeed perturbative small in weak coupling. 

Furthermore the vector bosons themselves acquire self energy corrections through fermion loops and vector boson 
loops through the non-abelian self- interaction. These are displayed in fig. H13(l . 




■^AAAy 




a 




^AAAi /WW 



FIG. 13: One loop diagrams contributing to the self-energy of charged and neutral vector bosons. 



The self-energy diagrams with fermion loops (diagrams (c) and (d) as well as similar diagrams with quarks for 
the neutral vector boson) yield a typical (HTL) self-energy correction of order 5^T^j31„^J and are perturbative for 
external momentum q ^ gT or for soft external momentum if M{T) ^ gT, where M{T) stands generally for the 
vector boson mass. 

These self energy corrections yield a Debye (electric) mass tod ~ oT [3iII3^I3^ . In non-abelian theories, a magnetic 
mass mm ^ g^T is also generated at high temperature fsil Is^ . 
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Since we are considering the high temperature hmit with T M\y,z{T) these contributions are perturbative for 
soft external momentum in the regime 



T:^ M{T):^ gT (V.l) 

A tighter bound emerges from consideration of the self-energy contribution from non-abelian loops, such as (a) 
and (b) in fig. H13|) . The (HTL) power counting [Sll l33 | for these diagrams is modified by the large momentum 
behavior of the propagators in the (physical) unitary gauge. Since for large momentum the vector boson propagators 
~ qi^q^/q'^]VP{T) ^ 1/M^(T) the naive (HTL) counting assigns a factor of 1/A'P{T) for each vector boson 
line, rather than the usual l/T^jsillHl. 

Therefore, according to this modified (HTL) power counting the vector boson loop with cubic (momentum depen- 
dent) vertex in diagram (a) in fig. H13|l. would yield a self energy contribution 



Four powers of temperature from the loop with two extra powers from the momentum dependence of each vertex 
and four powers of M in the denominator because each vector boson line yields a factor 1/M^(T). Similarly the 
vector boson loop with quartic (momentum independent) vertex in diagram (b) in fig. (|13|l . would yield a self energy 
contribution 



It is clear that if the usual (HTL) power counting holds for these diagrams as M{T) the infrared divergences 
lead to a breakdown of the perturbative expansion, requiring a resummation. 

However for sufficiently high temperatures but below the symmetry restoration scale Tew^ because the standard 
model features a smooth crossover from the broken to the unbroken "phase" , there is a regime for which Tew > T > 
M{T) within which a perturbative treatment is reliable. 

Within this window of validity of perturbation theory, the validity of the (HTL) approximation will be further bound 
by the region in which the self-energy contribution to the gauge vector bosons is perturbatively small. The analysis 
above reveals that for T ^ M{T) diagram (a) in fig. (|13|l gives the largest contribution, since at high temperature the 
loop momentum is g ~ T the self energy contribution HV.2fl will be perturbative if 



T > j^j^ =^ Af (T) > g^T (V.4) 

A tighter and more conservative bound is obtained from the observation that in using the free field theory propaga- 
tors for the internal lines the dispersion relation of on-shell vector bosons has been used. In particular the polarization 
vectors involve the bare mass ~ M{T). Hence a tighter bound is obtained by requiring 



^^^^^ ^ 'Wht) ^ ^'^^^^ ^ ^^-^^ 

The bounds above are the result of assuming the validity of the naive argument that combines the (HTL) power 
counting with the large momentum behavior of the vector boson propagator. 

However this simple counting argument ignores the subtleties associated with the underlying gauge invariance as 
well as potential cancellations. It is plausible that the term of 0{g^T^ /M'^{T) vanishes identically and that there is a 
cancellation of the terms of 0{g'^T'^ /]\-P{T) between the diagrams (a) and (b) as a result of Ward identities. // such 
is the case, then only a contribution of 0{g^T^) would remain in which case the perturbative expansion is reliable for 



M{T) > gT (V.6) 

An indication that these cancellations are plausible is found precisely in the neutrino self-energy in the high tem- 
perature limit obtained in the appendix and section (|III|I . The naive (HTL) power counting combined with the 
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large momentum behavior of the vector boson propagator cx \/M'^{T) would lead to the estimate ^ g '^T^/AP{T) 
for the first two terms and g'^T'^/M^k for the last two terms in the self-energy corrections (Ilii.llllli.l2|l ). However, 
remarkable cancellations in both terms lead to g'^T^uj/M'^(T) for the first two terms and g^T^/k for the last two 
terms, and perturbation theory is reliable for Af^(T) ^ gT. Of course, whether these cancellations also occur in the 
self-energy of the vector bosons must be studied in detail. 

In either of these cases (jV.4IV.6ll the vector boson masses from symmetry breaking M{T) are much greater than 
the Debye or magnetic masses, rriD ^ gT, mm ^ g^T respectively. 

Even if the more conservative bound (jV.Sp survives deeper scrutiny, within strict perturbation theory a regime in 
which the (HTL) approximation is reliable within this bound clearly exists, namely T > M{T) > g^T . However, this 
analysis warrants a deeper study of the screening corrections to the vector bosons to provide a reliable bound for the 
validity of the perturbative expansion at high temperature. 

While the (HTL) approximation is well understood |3lL l3^ in the absence of symmetry breaking, the interplay 
between the generation of a Debye and magnetic masses and the mass from symmetry breaking M(T) below the 
critical temperature is not yet completely understood. 

In ref.|33 the Debye screening corrections to the vector boson masses was studied in the temperature regime 
M{T) <C T <C VT2 w(0), where v(0) is the expectation value of the neutral component of the Higgs doublet at zero 
temperature. At finite temperature it is found in this reference that 

Therefore the regime of validity of the (HTL) approximation as stated in ref . [s^ corresponds to temperatures much 
smaller than the critical. 

In this regime of temperature it was found in ref. |33| that the vector bosons feature a longitudinal squared mass of 
the form 

Mf{T) ^ M^(T) + g^T^A{e^) (V.8) 
where A{9w) is a simple function of the Weinberg angle, and a transverse mass squared 

M^{T) = M^{T) (V.9) 

A magnetic mass term is expected at order (7^0,1131. These results are valid well below the critical temperature, 
and to the best of our knowledge there is as yet no clear understanding of the validity of the (HTL) approximation 
near the critical temperature. Ref. |33^ concludes by stating that a resummation program must be implemented, a 
statement that becomes more relevant near the critical temperature. This discussion is meant to bring to the fore the 
relevant but largely unexplored question of the validity of the (HTL) near the critical temperature. 

Such a study is obviously beyond the realm of this article and is deferred to future investigations. 

VI. DISCUSSION OF THE RESULTS 

• Collective modes: There are several remarkable differences between the spectrum of collective modes obtained 
above and those for fermionic collective excitations in QCD or QEDp^ IT^ Isil l33 |. To begin with, in the 
region in which perturbation theory is reliable, T ^ M\y,z{T) ^ gT, the spectrum depends on the mass 
parameter = gT/4, which determines the chirally symmetric gaps in the spectra, as well as the dimensionless 
ratio A = M^{T) /2m^T and features in general several branches. Gapped branches for quasiparticles and 
quasiholes are present for A < Ac ~ 1.275- •• a gapless branch for quasiholes exists for A < 7r/2 which 
becomes a gapless branch of quasiparticles for A > 7r/2. The quasihole branches terminate at end points that 
depend on the value of A and for large momentum k ^ TOj^A there are no quasihole branches available. For 
Ac < A < 7r/2 the quasiparticle spectrum features a pitchfork bifurcation with two branches emerging, the one 
with decreasing frequency terminates at an end point but the other with increasing frequency continues and 
merges asymptotically with the free field dispersion relation. For A > 7r/2 the collective modes are gapless 
quasiparticles whose dispersion relation lies below the light cone for k < vriv and approaches the free field 
dispersion relation for k 3> m^. 
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• Gauge invariance: An important aspect that must be discussed is the issue of gauge invariance. We have 
obtained the neutrino self-energy in the unitary gauge. This is a physical gauge in the sense that it displays only 
the physical excitations. In any other covariant gauge there are unphysical deg rees of freedom, in particular 
unphysical Goldstone bosons with a Yukawa coupling to neutrinos. In ref. 25] the neutrino self energy (at 
temperatures much smaller than the mass of the vector boson) was obtained in general covariant gauges an it is 
shown explicitly that the dispersion relation is independent of the gauge parameter. Furthermore, as discussed 
explicitly in references |17l l3lL 13^ the HTL approximation is gauge invariant. While this body of work clearly 
points out that the results obtained here are gauge invariant, there is a rather simple argument that makes the 
gauge invariance manifest up to the one loop order considered here. In the unitary gauge the only singularities 
in the Feynman propagators correspond to the on-shell propagation of physical degrees of freedom and the 
unitarity of the S-matrix follows directly from the Cutkosky rules 35] . In particular the imaginary part of the 
one-loop self-energy can be obtained directly from the Cutkosky rules at tree level, and these only involve the 
propagation of physical degrees of freedom. Since the full one-loop self-energy is obtained from a dispersion 
relation, the gauge invariance of the imaginary part guarantees the gauge invariance of the self-energy up to the 
order considered. For example, the calculation of decay rates in the standard model in the Born approximation 
in the unitary gauge yield the physical result 35. 41]. and the imaginary part of one- loop diagrams is obtained 
from the Born transition elements by the unitarity of the S-matrix, which is manifest in the unitary gauge [ssj 
(in fact this is the bonus and defining property of the unitary gauge). Therefore we conclude that the results 
obtained here are manifestly gauge invariant. 

• Width of collective excitations from vector boson decay: Another remarkable aspect that must be 
highlighted is that the width of the neutrino excitations in the medium is a result of the decay of the vector 
bosons into neutrino or neutrino- lepton pairs. This novel mechanism is rather different from the usual coUisional 
relaxation that leads to a width of the quasiparticles. In particular coUisional relaxation results in a width of 
order G\ since the Born amplitude for such process must necessarily involve the exchange of a vector boson 
and is therefore is of order Gp- In contrast, the width acquired via the decay of a more massive vector boson 
is of order Gp since the Born amplitude is of order \/Gf- We emphasize that this mechanism is present even if 
the vector boson population in the medium is negligible, as can be gleaned from the fact that the term from the 
imaginary part of the self energy with support on the mass shell of the neutrino excitation involves Nf + Nf,. 
Therefore even if iV^ there will be a width for the neutrinos if there are fermions (either neutrinos or leptons) 
in the medium. There is a simple interpretation of this rather striking result, which can be gleaned from the 
discussion of the width in terms of the kinetic equation ljIV.40|) . While the first term in (jIV.40|l which describes 
the gain in neutrino population from the decay of the vector boson, vanishes when Nyy ~ 0, the second term 
does not vanish when Nyy ^ 0. This second term describes the recombination process in which the neutrino 
in the medium annihilates with either another neutrino (neutral current) or a lepton (charged current) in the 
medium to produce a vector boson, the stimulated emission makes this term non-vanishing even for vanishing 
population of vector bosons. Therefore this contribution to the width of neutrino excitations will always be 
available in the medium provided the kinematics for energy-momentum conservation is fulfilled. Indeed this is 
the major restriction for the width that arises from this process, in order for the recombination process to be 
available, energy-momentum conservation require that the momentum of the lepton (or neutrino) in the bath 
must be very large (when the mass of the vector boson is much larger than that of the neutrino quasiparticle or 
quasihole) and therefore this process will probe the tail of the fermionic distribution function with the ensuing 
exponential suppression at low temperatures. This novel relaxational mechanism has been previously studied 
in different contexts in refs. [3^l40| . 

• Comparison with coUisional damping: In references pol |2^ the coUisional damping rate of neutrinos was 
computed up to two loops order for temperatures much smaller than the vector boson mass. While we cannot 
directly compare those results to ours for two reasons: a) our study focuses on temperatures much larger than 
the (temperature dependent) vector boson mass, b) we have obtained the damping rate of quasiparticles and 
quasiholes rather than weakly interacting neutrino states, we can compare the results to obtain at least an 
estimate of the new effects. The result for the coUisional damping rate obtained in refs.|23,|23 is 



re(w) - GlT^LU (VI.l) 

where lu is the energy of the neutrino, uj ^ k for a fast moving neutrino in the medium. In ref. 27j a similar result 
for the coUisional width was obtained but temperature T replaces the neutrino energy in eqn. (|VI.1|) . These 
results are obtained from a kinetic equation that includes the different scattering contributions, or alternatively 
as in ref. ii2,(|] from an analysis of the imaginary part of the neutrino self-energy up to two loops. These calculations 
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take the external neutrino to be described by a free field state and are restricted to temperatures below the 
vector boson mass. Since fast moving quasiparticles are very similar to the weakly interacting neutrino states, we 
can proceed to compare the results for the width of the quasiparticle excitations for large momentum k ^ A 
given by eqn. (|IV.56|I with the result of eqn. (|VI.ip for lo ^ k. The ratio of these results is 



r+(fc) r 4fc 1 

r,(fc) ^ Ang^ V k \ Lm^AJ ^ ' 



In the expression above we have (unjustifiably) taken the vector boson mass in (jVI.ip to be equal to M{T) in 
our calculation, with the purpose of offering a comparison. Granting these caveats, it is clear that under the 
perturbative assumption M{T) ^ m,, ~ gT the ratio r+(fc)/rc(fc) > 1 for 1 ^ 2i|A > gs. Hence, for a fairly 
wide range of neutrino energies, the relaxation mechanism via the decay of vector bosons could be comparable 
to that via collisions. While this comparison must be considered only an estimate (given the caveats mentioned 
above), it at least suggests that the damping mechanism found in this article is competitive with if not larger than 
the usual coUisional relaxation and must be included in treatments of non-equilibrium phenomena of neutrinos 
in the early Universe. 

• Regime of validity: Although the results of ref.jH3| are valid well below the critical temperature, they are 
encouraging in that the vector bosons acquire a Debye screening mass of order gT which, if taken at face value 
near the critical temperature would entail that A ^ g 1. In this case the dispersion relation of the collective 
modes will feature the wealth of novel phenomena studied above. If further studies confirm that perturbation 
theory is reliable very near the critical temperature and the vector boson acquires a screening mass of order 
gT then this would translate in that A ^ g. Therefore, at least within perturbation theory there is a wide 
range I > A > g for which the collective modes feature the most striking and richer aspects of the dispersion 
relations found above. A similar conclusion is reached if it is required that Af^(T) > g^T^ so that the transverse 
mass is larger than the magnetic mass to avoid the breakdown of perturbation theory at long wavelengths. 
If, on the other hand, the more stringent constraints l|V.4IV.5() are shown to be the relevant ones because of 
infrared phenomena associated with the magnetic mass very near the critical temperature, then A > 1 and 
the properties of the collective modes depart perturbatively from those of the vacuum states. It is clear that 
a further assessment requires a deeper understanding of the validity of the (HTL) program near the critical 
temperature. 



VII. CONCLUSIONS 



Motivated to explore non-equilibrium properties of neutrinos that could impact on thermal leptogenesis, or more 
generally in neutrino transport in the early Universe, we studied the collective excitations of neutrinos in the standard 
model at high temperature but below the symmetry breaking scale Tew- The main assumption is that the transition 
from the broken to the unbroken symmetry states in the standard model is either a smooth crossover, as supported 
by the lattice data with the current bound on the Higgs mass, or of second order. In this scenario the expectation 
value of the neutral scalar in the standard model vanishes continuously near the transition resulting in that the mass 
of the vector bosons becomes smaller, which in turn leads to a large population of vector bosons in the thermal bath. 
We have obtained the spectrum of collective excitations for standard model neutrinos in the (HTL) approximation 
in the regime T 3> Mw,ziT) S> gT within which perturbation theory is reliable. The excitation spectrum consists 
of left handed positive energy negative helicity quasiparticles and left handed positive energy and positive helicity 
quasiholes. Antiquasiparticles and antiquasiholes carry negative energy and opposite helicity assignments. The 
excitation spectrum features a chirally non-breaking mass scale 



- ^. (VII.l) 

and depends on the dimensionless ratio A = M^{T)/2m^T. For A < Ac = 1.275- ••, the quasiparticle spectrum 
features two gapped branches and the quasi-hole spectrum features one gapless and two gapped branches of collective 
modes. 

The lower quasiparticle branch as well as all the quasihole branches terminate at particular end point values of 
the momentum that depend on A. For A > 7r/2 only a gapless quasiparticle branch is available, with a dispersion 
relation that is below the light cone for k <^ mi, and asymptotically reaches the free particle dispersion relation for 
k > m^. 
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A novel feature revealed by this study is that the decay of the vector bosons into neutrinos and leptons leads to 
a damping rate for the neutrino collective modes. A simple kinetic interpretation of this phenomenon was given and 
analytic expressions for the damping rates were obtained in the limit of small and large momentum of the excitations. 
These are given by equations l|IV.48|) and (|IV.56|I respectively. For A <C 1 quasiparticles and quasiholes at rest in the 
lower branch are heavily damped while those in the upper branch are weakly damped, while for A > 1 all collective 
modes at rest are heavily damped. Fast moving quasiparticles are weakly damped. 

We have compared the damping rate of fast moving quasiparticles to the results for the coUisional damping rate of 
refs.|20, 26, 27] which is a two loop result {0{Gjp)), the ratio is given by equation (jVI.2p . The damping rate resulting 
from vector boson decay at one loop can be larger than the coUisional relaxation rate in a wide range of neutrino 
energy. 

We have discussed in detail the gauge invariance of the results and assessed the region of validity of the (HTL) 
approximation which relies on a strict perturbative expansion. 

The novel aspects of the neutrino collective modes near the crossover (or second order) transition to the unbroken 
phase as well as their non-equilibrium relaxational dynamics studied in this article could be important for a reliable 
assessment of mechanisms for thermal leptogenesis 

Further questions: Clearly an aspect that warrants further study is a more quantitative assessment of the 
screening corrections to the gauge boson propagators, namely the diagrams displayed in fig. in particular 

diagrams (a,b). A detailed evaluation of these self-energy corrections will yield a better assessment of the domain 
of validity of the (HTL) approximation near the critical temperature. At temperatures closer to the transition (or 
crossover) a resummation of the perturbative series must be invoked in order to correctly treat the infrared divergences 
in the vector boson propagatorsj42i| . 

While we have focused on temperatures below the electroweak scale, the study of neutrino collective excitations 
above this scale is also inherently important for a deeper understanding of early Universe neutrino cosmology. Above 
this scale, in the phase where the SU{2) (g) U{\) is restored, new interaction vertices between neutrino, leptons and 
the neutral scalar emerge. The contribution from these vertices to the collective excitations and their damping rates 
is an important future direction of study. 

This analysis will require a deeper understanding of the validity of the perturbative expansion as well as the issue 
of gauge invariance as higher order corrections are contemplated. 

Finally, going beyond the standard model, it becomes important to study the effect of these novel phenomena 
upon neutrino oscillation and in particular the possibility of MSW resonances. An intriguing possibility is that of 
oscillations between different branches as well as the competition between the oscillation and relaxation time scales. 
Furthermore small neutrino masses will lead to helicity flip transitions which may in turn lead to radiative transitions 
between the quasiparticle and quasihole branches (or their antiparticles). 

These are currently the focus of our current studies on which we expect to report in a forthcoming article ,34J. 
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APPENDIX A: REAL TIME PROPAGATORS AND SELF ENERGIES. 

1 . Fermions 

For a generic massless fermion field /(x, t) the Wightmann and Green's functions at finite temperature are the 
following 

iSl^{x-x',t-t') = (/„(f,t)7^(f' ,0) ^ ^Y."^^^'^^'"'^ ^Slp{P^t-t') (A.l) 

p 

is<p{x - f , t - 1') = -(7^(x', t')Ux;, t))^^Y. ^''^'^^ '^^Ap^ ^ - 

p 

where a, (3 are Dirac indices and V is the quantization volume. 
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The real time Green's functions along the forward (+) and backward (— ) branches are given in terms of these 
Wightmann functions as follows 

{f^+\x,t)f^\x',t')) = iS++{x-x',t-t')=iS>{x-x',t-f')Q{t-t') + iS<{x-x',t-t')Q{t' -t) (A.3) 
{f^+\x,t)f^\x' ,f')) - iS+-{x-x-',t-t') = iS<{x-x',t-t') (A.4) 

At finite temperature T it is straightforward to obtain these correlation functions by expanding the free Fermion 
fields in terms of Fock creation and annihilation operators and massless spinors. The result is conveniently written in 
a dispersive form 

/oo 
^ipoP>;3(p,Po)e-*^°(*-*') (A.5) 
-OO 

/OO 
dpnpl0{p,Po)e-'^''^'-''^ (A.6) 
-OO 

where 

P^Ap^Po) = [l-^/bo)]pl/3(P.Po) (A.7) 
PaAP^Po) = ^fiPo)pi,o{P,Po)- (A.8) 
and Nf{po) is the Fermi-Dirac distribution function 

Nfipo) = (A.9) 
e T + 1 

The free Fermion spectral density is given by 

P^{p,Po) = ^Sipo-p) + ^5{po+p) (A.IO) 
li± = 1°PT1-P (A.ll) 

2. Vector bosons 

Consider a generic real vector boson field An{x,t) of mass M. In unitary gauge it can be expanded in terms of 
Fock creation and annihilation operators of physical states with three polarizations as 

1 f'^(k) r - J. - 1 

= ^ V V^L^ ag;^e-*'"'=*e*'=-^ + at e*'"'=*e-*^ ; k^'ci.^xik) = (A.12) 



where cok = + and k^^ is on-shell k^ = [wk , k) . The three polarization vectors are such that 

J^^ikVxik) = P'^ik) = -U''-^) (A.13) 

A=l ^ ^ 

It is now straightforward to compute the Wightmann functions of the vector bosons in a state in which the physical 
degrees of freedom are in thermal equilibrium at temperature T. These are given by 

{A^{x,t)A,{x\t')) = tG>Jx-x',t-t') (A.14) 

{A,{x',t')A^{x,t)) = iG<ix-x',t-t') (A.15) 
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where G ■ can be conveniently written as spectral integrals in the form 



iGl„{x- x' ,t-t') 
iG< „{x- x',t-t') 



V 



V 



(A.16) 
(A.17) 



where 



mko) = — 

e T — 1 



and the spectral density is given by 



P^iv{ki3, k) = - — P^ii,(k)5(kn - LOk) - P^ii,{-k) d{ko + uJk) 

ZUJi, L 



(A.18) 



(A.19) 



In terms of these Wightmann functions the real time correlation functions along the forward and backward time 
branches are given by 



(4+)(f,t)4+)(x',t')) = iGl,{x-x\t-t')Q{t-t')+zGl,{x-x\t-t')Q{t' -t) (A.20) 
(4+)(f,i)4-)(x',t')) = iG<{x-x;',t-t') (A.21) 



For the charged vector bosons the correlation functions can be found simply from those of the real vector boson 
fields described above by writing the charged fields as linear combinations of two real fields A^^^, namely 



W^{x-,t) ^ ^{Al{x,t)±iAl{x,t)) 



(A.22) 



It is straightforward to find the correlation function 



{W+{x, t)W- {x', t')) - G>,{x -x\t~ t') 
and similarly for the other necessary Wightmann and Green's functions. 



(A.23) 



3. Retarded self-energies 

The diagrams for the one-loop retarded self-energy from charged current interactions are displayed in fig. H14|l . A 
straightforward calculation yields for the charged current contribution the following result 



S^^?(f -x',t- t') = -Y^l^ [iS++{x -x',t^ t')iG++{x -x',t- t') - iS<{x -x',t- t')iG<^{x -x',t- t')] Y L 

(A.24) 



with 



(1 + 7^ . (1^7^) 
2 ' 2 ■ 

A similar result is obtained for the neutral current contribution to the self energy by simply replacing 
g/2cos9w and Mw Mz = Mw/ cos 9^. 
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x't') 




FIG. 14: Retarded self-energy for charged current interactions. The wiggly hne is a charged vector boson and the dashed line 
a lepton. The labels (±) correspond to the forward (+) and backward ( — ) time branches. The corresponding propagators are 
iS^'^ {x — x',t — t') and iG^^{x — x',t — t') for leptons and charged bosons respectively. 



Using the representation of the fermion and vector boson propagators given above the retarded self-energy ljA.24|l 
can be written as 



T,ret{x -x\t^t') ^ I R \^w{k, h) + T.z{k, fco)l L e»'^-(^-=^') e-''=«(*-*')e(i - t'), (A.25) 



The contribution from charged and neutral vector bosons are given by 



T.w{k, fco) 



.2 , ^3^ 



2 J (27r)3 



dpo / dqo5{p^ + qo- ko)-f^pf{k-q,po)p%ik)iq,qo)Y (1 - Nfipo) + Ntiqo)) (A.26) 



dpo / dqoS{po + qo~ko)-f''pfik-q,po)pf,M<lo)l''{'i^-Nfipa) + Nhiqo)) (A.27) 



4cos"^ Ow J {2T:y 

where pw,ziq, qo) are the vector boson spectral densities given by IjA.lQp with M = Mw,z{T) respectively. It is clear 
that Yiw,z{k, fco)) correspond to a vector-like theory. 

Using the integral representation of the function Q(t ~ t') the retarded self energy can be written in the following 
simple dispersive form 



T.ret{x - x' ,t - = 

k 



2tt 



Y.w{k,uj) + Y.z{k,uj] 



L 



(A.28) 



I]n/,z(fc,w) = / dfco 



^w,z{k, fcp) 
kt) — LO — ie 



(A.29) 



where e — > 0"*". 

Hence from the above expression we identify 



Ysw,z{k,uj) = -ImY.w.z{k,u!) 



(A.30) 



Furthermore, since we are considering massless neutrinos and leptons in the high temperature approximation, the 
fermionic spectral density is proportional to 7 matrices and does not feature the identity matrix or 7^ (see eqns. 
(|A.10p . (jA.lip therefore there is the following simplification 



R Y,w{k,uj) + ^z{k,oj) L~ J^w{k,uj) + Ilz{k,oj) L 



(A.31) 
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